Phonon Coherence and New Set of Sidebands in Phonon- Assisted Photoluminescence 
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We investigate excitonic polaron states comprising a local exciton and phonons in the longitudinal 
optical (LO) mode by solving the Schrodinger equation. We derive an exact expression for the 
ground state (GS), which includes multi-phonon components with coefficients satisfying the Huang- 
Rhys factors. The recombination of GS and excited polaron states gives one set of sidebands in 
photoluminescence (PL): the multi-phonon components in the GS produce the Stokes lines and 
the zero-phonon components in the excited states produce the anti-Stokes lines. By introducing 
the mixing of the LO mode and environal phonon modes, the exciton will also couple with the 
latter, and the resultant polaron states result in another set of phonon sidebands. This set has a 
zero-phonon line higher and wider than that of the first set due to the tremendous number of the 
environal modes. The energy spacing between the zero-phonon lines of the first and second sets is 
proved to be the binding energy of the GS state. The common exciton origin of these two sets can 
be further verified by a characteristic Fano lineshape induced by the coherence in the mixing of the 
LO and the environal modes. 

PACS numbers: 78.20.Bh; 42.50.Ct; 78.55.Et 
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It has been recognized that the interactions of ex- 
citons with phonons in light-emitting materials lead 
to a number of interesting effects in optical proper- 
ties, such as the phonon sidebands in photoluminescence 

SL) spectra, the exciton dephasing and self-trapping 
S II S H H H H . Theoretically, the sideband struc- 
tures caused by electron-phonon interactions for local ex- 
citons at deep centers were investigated using the Huang- 
Rhys model 0, 0, 01 . In these theories, an adiabatic 
approximation is applied to separate the degrees of free- 
dom of excitons and phonons. One set of Stokes lines 
(SLs) in PL, observed in the experiments, has been very 
successfully explained by the theories. However, there are 
some features in the experimental spectra associated with 
the interactions of excitons and phonons, such as compli- 
cated structure containing more than one set of phonon 
sidebands, anti-Stokes lines (ASLs), and Fano-like line- 
shape, seemingly beyond the adiabatic approximation. 
Recently, the fabrication of nanostructures has renewed 
the study of this field, as the confinement in quantum 
dots much enhances the interaction between carriers and 
phonons. Such an enhancement has been experimentally 
observed jljj, |13| . Fomin et al. have proposed a nona- 
diabatic approach to interpret the unusual phonon side- 
bands observed in quantum dots Verzelen et al. 
have indicated that the exciton-phonon bound states as 
a whole, called excitonic polarons, and the decaying of 
their phonon components into phonon thermostat have 
to be put under consideration in cases of strong exciton- 
phonon interactions [l5| . 

In this Letter, we present a unified theory beyond the 
adiabatic approximation to explain complicated struc- 
tures of phonon sidebands in PL. The exact eigenstates 
and eigenenergies of excitonic polarons are obtained by 
solving the Schrodinger equation of an exciton interact- 



ing with LO mode. The obtained ground state (GS) con- 
tains multi-phonon components with coefficients strictly 
satisfying the Huang-Rhys factors, which gives the SLs 
including a ZPL. Meanwhile the zero-phonon components 
of the excited states can produce the ASLs. These SLs 
and ASLs form one set of sidebands. By introducing 
the mixing of LO mode and environal modes, another 
set emerges. The energy shift between two sets is deter- 
mined by the strength of exciton-LO phonon coupling, 
and the coherence in the mode mixing results in specific 
Fano-like lineshape. 

First, let us consider the interaction between an exci- 
ton and a LO mode 

H = eoa'^a + hcjoblbo + Voa''a{bl + bo), (1) 

where a and bo are annihilation operators for exciton and 
for LO phonon, eg and Tiujo are their energies, respec- 
tively, and Vo is the coupling strength between exciton 
and phonon. Here we regard the exciton as one quasi- 
particle, supposing that it has a large binding energy. 

For Hamiltonian , the mth eigen-wavefunctions can 
be written as 
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where n is the number of phonons of the LO mode in 
the corresponding component, |0) is the vacuum, and the 
coefficients satisfy the following iteration relations 



[E,n - ifO+nhuJo)]Crn:n = Voy/nC,n:n-l+VoVn + lCm.;n+l- 

(3) 

For the GS {m — 0), Eo is below eo, and coefficients are 
related with each other by 

(£^0 — eo)co;0 
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C0;2 = 



{Eq - £o - Jii^o)co-i - VqCq-o 
V2Vo 
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From these relations, we can find the exact expressions 
of the GS energy 



Eq — £0 — Aq 



(5) 



with Aq = j-^ being the binding energy of the excitonic 
polaron, and the coefficients 



CO;n 
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(6) 



T/2 A 

with 5 being the Huang-Rhys factor S — = 
Considering that the nth component in GS '^m=o is as- 
sociated with n phonons, |co;nP gives the intensity of a 
phonon sideband at energy Eq — nhujQ. This is as pre- 
dicted by the Huang-Rhys theory. 
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FIG. 1: Several low-lying eigenenergies Em as functions of 
'^max, the maximum number of phonons included in the cal- 
culation, eo = 3.366eV, hcoo = 0.069068eV, (a) Vo = O.OleV, 
(b) Vo = O.OSeV, (c) Vo = 0.05eV, (d) Vo = O.OSeV. The dot- 
ted lines represent the corresponding values of eo + mh,u)o for 
comparison. 

For excited states '^m with m > 0, exact analytical 
expressions can not be obtained in this way. However, 



from a numerical diagonalization shown in Fig. 1 one 
can see that for both cases of weak and strong coupling 
the eigenvalues can be expressed as 



Em = mflLLlo - Aq. 



(7) 



From this the coefficients in the eigenfunctions can be 
obtained by the foUowing iteration relations 
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where 
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for n = 0,1,2, 
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Aq — (m — n — l)?icLio — [n 
with the initial value 

1 
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Aq — mfiujo ' 

and by the normalization. Although the eigenenergy of 
state with m > is higher than Eq, it should have a 
finite occupation rate as long as the photon energy of the 
excitation laser is high enough . In order to calculate 
the PL spectrum of the polaron states, we must know 
the population of each state. For simplicity, we assume 
that all the states with energies less than the laser energy 
^-lascr have the same occupation rate. The intensity of PL 
at photon energy e is proportional to 



^6l(£:iascr - Em)S{e ~ Em + uflUJo). 



(9) 

The luminescence peaks produced by the GS ^0 con- 
sist of a ZPL at energy Eq = eq — Aq and a number of 
SLs at energies Eq — nhujQ. This ZPL originates from the 
zero-phonon component (cq-q) of the GS state while the 
remaining lines are from the components {co-n>a) with 
more than one phonon. The intensities of the SLs rela- 
tive to the ZPL follow a Poisson distribution predicted 
by the Huang-Rhys theory. Other eigenstates ^'m>o will 
produce complex sideband structures, which is not dealt 
with in the Huang-Rhys theory. Radiative recombination 
of these states will intensify the ZPL peak by accompany- 
ing with simultaneous emission of m LO phonons. When 
they emit n > m phonons in the radiative recombination, 
the intensity of the SL at — nfiojQ will be enhanced. 
Meanwhile, the component Cm-n<m of state may pro- 
duce the ASL peak at Em—nhcuQ above the ZPL, together 
with the emission of n(< m) phonons. The relative in- 
tensity of this ASL is Pm;m-nl Pm-TU = Ji{^-n)\ - 

It is known that besides LO mode, there are a large 
number of bath phonon modes in solids. Although these 
bath modes do not directly interact with the exciton, 
they may have crucial influence on phonon components 
of the polaron states. This influence can be accounted 
for using the following Hamiltonian of the mode mixing 

^1 = E 5(&S^A + H.c.) + hioxiblbx + i), (10) 

X X 
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where A is the mode index of bath phonons, and g is the 
strength of the mixing. For simpUcity, g is assumed to be 
the same for all bath modes and the dispersion relation 
ujx is continuous with a constant density of states between 
energies and W. For a pure harmonic phonon system 
without exciton, the mixing terms should vanish as all 
the phonon modes are orthogonal to each other. The 
formation of the excitonic polaron states, however, vio- 
lates the orthogonality between the LO and bath modes, 
and thus creates the mixing terms. It is this mixing that 
makes the LO phonons in the polaron states decays 
into the continuum of bath modes. From Hamiltonian 
Hl()|l . the orthogonal modes become 



huJi - hujo 
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with frequency coi determined by solving the equation 

^2 



huJi = hujQ + ^ 
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(11) 

m 
(12) 



The resultant set of orthogonal modes {/} will couples 
with the exciton at a strength 



Vi 



(13) 



Figure 2 shows the phonon frequency dependence of Vi . 
It can be seen that most of the modes {/} are still weakly 
coupled to the exciton except for those modes with fre- 
quencies LOi ~ ojQ. As mentioned earlier, this results in 
not only the broadening of the peaks in the first set but 
also the creation of the second set of the phonon side- 
bands. 
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strongly couple with the exciton; P2 includes the re- 
maining modes weakly coupled to the exciton. Because 
the modes in Pi with frequencies in the vicinity of wq 
have coupling strength of about Vb with the exciton, the 
formed polaron states will produce the first set of phonon 
sidebands including broadened ZPL, SLs and ASLs, as al- 
ready discussed above. It should be noted that the ZPL 
is centered at the GS eigenenergy Eq which lows with 
increasing g due to an increase of the mode number in 
group Pi as seen from Fig. 2. Now we turn to discuss the 
modes in P2. For a mode / in this group, it has a weak 
coupling with the exciton, i.e., Vi 0, so the binding 
energy of the relevant polaron state is approaching zero. 
In other words, the eigenenergy of the polaron state is 
very close to the bare exciton energy cq. Thus, the ZPL 
produced by the zero-phonon component in this state al- 
most located at eo, higher than the ZPL of the first set by 
an energy about Eq — Eq. The intensity of the ZPL in the 
second set should be much stronger than that of the ZPL 
in the first set since the ZPL intensity is proportional to 
the number of involved phonon modes. Because of the 
evident dispersion of the bath modes in group P2, it can 
be expected that the ZPL of the second set is wider than 
that of the first set. The energetic position of the one- 
phonon sideband in the second set ranges from eo (the 
ZPL position) to eo — W. The transition probability of 
this sideband is found to be 

{huJi — huJo + g)^ 



g2 + (fiujQ _ fioJiY 



(14) 



From Eq. (|14|1 . it can be easily known that pm.i reaches 
its maximum value for huJi — Hujq + g and the minimum 
for tiuji — hujQ — g. When the phonon frequency is in- 
creased from huo — g to huJo + g the intensity of the 
sideband rapidly increases. However, as the phonon fre- 
quency is further increased over Hujq + g the sideband 
intensity slowly decreases. This results in a typical Fano 
asymmetric lineshape. As a matter of fact, one can 
rewrite Eq. 114|) into a standard Fano lineshape func- 
tion a; = (e -I- l)^/(e^ + 1) with e = {Hloi — hiOnUg. Here 
the Fano asymmetry parameter is the unity jl7|. Since 
the mixing between LO-mode and A-mode (which forms 
the I mode components in the polaron states, illustrated 
in Eq. (|ll|l l. introduced by -ffi, is coherent, an interfer- 
ence between the optical transitions associated with the 
two phonon components in the / mode will occur. That 
is the origin of Fano interference in the phonon-assisted 
luminescence. 

By diagonalizing the whole Hamiltonian H + Hi , one 
can calculate the PL spectrum by 



FIG. 2: Coupling strength Vi as a function of phonon fre- p(e) 
quency. 
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To give a clarification, we can divide the modes {/} 
into two groups: Pi includes modes with uji ^ luq which 



X 5{e — E„i + riifiuji 



(15) 
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FIG. 3: Relative intensity as a function of photon energy 
when Hi is included. In the simulation 1000 modes uniformly 
distributed in a band from to O.leV are included. Vo = 
0.0078eV. Other parameters are the same as those in Fig. 1. 

where Em and Cm-{ni} are eigenenergy and eigenfunction 
of the whole Hamiltonian. In Fig. 3, we plot the PL 
spectrum obtained from a numerical calculation. As ex- 
pected, there are two ZPL peaks, 0' and 0, belonging to 
the first and second set, respectively. In each set, there 
are several SLs and ASLs. The SLs are red-shifted from 
their corresponding ZPL by integer times Tiwo while the 
energy distances of the ASLs from the ZPL are not inte- 
ger times hujQ. It is because that the energy positions of 
the ASLs are determined by energies of the excited po- 
laron states whose energies are not exactly higher than 
the GS energy by integer times Tiloq. The Fano interfer- 
ence features are partially smeared by the summation in 
Eq. H15|l . but they can still be seen near some peaks. 
These results catch the essential nature observed in ex- 
periments [lij . 

In a summary, by reconsidering the exciton-phonon in- 
teraction and decaying of LO phonons we can conclude 
that there are two sets of phonon assisted peaks: one is 
from excitonic polaron states of the phonon modes with 



frequencies loi ~ wq, and the other is from polaron states 
of the remaining bath phonon modes. The latter set was 
not properly considered in previous theories and was usu- 
ally assigned to be the first set because it may include 
the strongest "ZPL" peak. The common exciton origin 
of these two sets of phonon sidebands can be further veri- 
fied by the characteristic energy difference between them 
determined by the exciton-LO phonon coupling strength, 
and by the specific Fano lineshape reflecting the coher- 
ence in the phonon mode mixing. Together with the ex- 
perimental observation, the present theory may give a 
new and better understanding of the complicated struc- 
tures in the phonon-assisted luminescence spectrum. 
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